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Imaging  through  Turbulence 


Typical  Imaging  Conditions  at  0.5  |im  at  AMOS 


Simulations  of  the  Hubble  Space  Telescope  as  it  would 
appear  from  the  3.6  m  AEOS  telescope  at  a  range  of  700 
km  in  1  ms  exposures  at  0.9  pm  wavelength  under  a 
range  of  seeing  conditions. 


Pristine 


Good  Seeing 
Terminator 


Marginal  Seeing 
Terminator 


Daytime 
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PROTEA 

<#•#  (Partnership  for  Research  in  Optical  Technology  and  Estimation-theoretic  Analysis) 


*  New  image  restoration  methods  to  provide  24-hour 
horizon-to-horizon  coverage 

*  Information  theoretic  approach  to  Space  Object 
Identification 

*  Pre-conditioning  of  data  to  reduce  computational 
overhead 

*  Understanding  speckle  phenomenology  in 
Superresolution  Active  Modulation  Imaging  (SAMI) 

*  Modeling  information  communicated  through  the 
image  chain 
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Algorithms  for  Multi-Frame 
Blind  Deconvolution 


Approaches  for  MFBD  algorithms  that: 

•  Enable  faster  convergence  with  higher  likelihood  of  achieving  the 
global  minimum. 

•  Improve  quality  of  image  restorations. 

•  Extend  limits  of  SSA  imaging  capabilities. 


D/rO  =  45 


Truth 


New  FF  alg  Standard  alg 


James  Nagy,  Emory  University 
Stuart  Jefferies,  U  Hawaii 


Initial  PSF  Estimates  using  Wave  Front 

Sensor 

•  WFS  collects  coarse  grid  gradients. 

•  Overlapping  measurements,  and  frozen 
flow  assumption  provides  finer  grid 
approximations. 

•  Novel  interpolation  and  constrained 
optimization  algorithms  reconstruct  initial 
estimate  of  PSF. 

•Substantial  improvement  over  adaptive 
optics  for  high  D/r0 
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Family  of  3-peak  Gaussian  Sources: 
Prior  on  Central  Peak  Location 


Uniform  prior  on  location 
of  central  peak 


Comparison  of  Quadratic  Regularizers  -  Suyu  Approach 


MSE  vs.  Regularizer  Strength  for  3-Peak  Gaussian  Scene 


Current  accommodation-based  metric  of 
optimality  -  picks  curvature  regularizer! 


Previous  approach  -  Evidence-based 
metric  also  picks  curvature  regularizer! 


-10 


-5 


log  x 
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Image  reconstruction  of  spacecraft  observed 
through  atmospheric  turbulence  is  limited  by 
unmodeled  effects  such  as  transient  glints: 


Consecutive  images  of  HST  at  10  Hz 


Robust  MFBD :  Algorithm  is  iteratively 
re  weighted  blind  deconvolution: 

1 .  Assign  weight  to  each  data  pixel 

2.  Reconstruct  image  using  weighted  data 

3.  Update  weights  where  predicted  image 
does  not  agree  with  data 

4.  Repeat. 


Transient  glints 


Simulated 
frame  with 
glint 


Traditional 

MFBD 


New  robust 

MFBD 


Unmodeled  noise 


Simulated  frame  Traditional  New  robust 
with  unmodeled  MFBD  MFBD 

noise 


Edge  effects 


V 


Simulated  frame  Traditional  New  robust 
clipped  at  edge  MFBD  MFBD 


Brandoch  Calef,  Boeing 
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High-resolution  imaging  of  LEO  targets 
through  stror^g  atmospheric  turbulence 


Raw  data  frame 


Single  telescope 
+MFBD 

+Spectral  nulls 


Multi-telescopes 

+Cross-channel 

consistency 


WFS  measurements 
+Frozen  Flow 
Model  (FFM) 


Multi-telescope 
+Annular  partition 
+Multi  Aperture 
Phase  Retrieval 


+  FFM  +  MFBD+WFS 


Increasina  svnerav  between  data  acauisition  and  orocessin 


Increasina  aperture  diversit 


Benefits  for  SSA:  extend  sky  coverage  and  range  of  observing  conditions 


Stuart  Jefferies,  University  of  Hawaii 
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Robust,  flexible  image  reconstruction: 
passive  sensing  of  objects  in  GEO 


•  Amplitude  interferometry 

•  Intensity  interferometry 

•  Shadow  imaging 


*  Quantify  performance  limits: 

•  SNR 

•  Cramer-Rao  Lower  Bounds 

•  Reconstructed  image  quality 

•  Etc. 


Amplitude 
Interferometry  (Al) 


Intensity 

Interferometry  (II) 


Shadow  Imaging 


PROS 


Higher  measurement  SNR 
compared  to  II,  potential  to 
measure  phase  information  in 
addition  to  Fourier  modulus 

Path  matching/stability 
requirements  less  stringent  than 
for  Al,  resolution  not  limited  by 
atmospheric  turbulence 

SNR  limited  by  background  star 
not  object  itself 


CONS 


Exacting  requirements  for  path  matching/stability, 
ground-based  version  impacted  by  atmospheric 
turbulence 

Fundamental  SNR  is  extremely  poor 


Measures  object  silhouette  only,  limited  number  of 
stars  means  limited  observation  opportunities, 
sensor  bandwidth  requirements  most  likely  increase 
going  to  a  space-based  implementation 
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Rotating  PSF  via  Spiraling  Pupil  Phase 


Sudhakar  Prasad,  University  of  New  Mexico 


Rotating  PSF  (for  M=7  Fresnel  zones) 


fo 

3 

o 

VO 

0/ 

+24  rad  +1 6  rad  +8  rad  0  rad  -8  rad  -1 6  rad  -24  rad 


Changing  Defocus  g 

+24  rad  +1 6  rad  +8  rad  0  rad  -8  rad  -1 6  rad  -24  rad 


Conventional  PSF  -  spreads  dramatically,  rapid  loss  of  sensitivity! 


•  Rotating  PSF  - 

•  Single-lobe  PSF  with  relatively  stable  shape/size 

•  High  3D  image  capture/reconstruction  sensitivity 

•  Defocus  g  =  Sz.Q/?i,  Q  =  R2/l02,  solid  angle  of  imager  at  source 

•  SSA  application  - 

•  Ranging  of  space  debris  relative  to  AF  assets 

•  Can  reconstruct  range  and  brightness  in  a  crowded  3D  scene 


Reduced-Measurement  Active  Polarimetry 


Full  Mueller  polarimetry  requires  at  least 
]  16  measurements. 

Reduced-measurement  polarimeter 
concepts  that  control  illumination  and 
analysis  polarization  states  to  maximize 
target  contrast 

(b)  Image  of  Spectrolon  from  SWIMMP 


Students  from  UA  constructed  an  active  1.5-jum  partial  polarimeter  at  WPAFB  for  use  in  a  basic  research  program 
on  human  signature  detection.  The  polarimeter  makes  at  most  12  measurements,  but  can  be  operated  with  a 
single  polarized  measurement.  Initial  lab  and  outdoor  testing  was  completed  in  the  fall  of  2012,  and  development 
will  continue  in  2013. 


Funded  in  2012  to  transition  concepts  from  our  AFOSR-sponsored  basic  reasearch  to: 

•  AFRL/RYJT  (signature  detection)  •  USA/NVESD  (IED  detection)  •  Sunkist  (Optical  Fruit  Grading) 
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Scott  Tyo,  U  Arizona 


Peregrine:  World's  First 
Collapsible  Space  Telescope 

*  Photon  Sieve  Observations  of  the  Sun 

-  Deployed  from  3U  CubeSat  as  part  of  FalconSAT-7 

-  0.2m  diameter,  f  =  0.4m  I  =  656. 5nm  (H-alpha),  2.5  billion  holes! 

-  Launch  date  of  early  2015,  total  lifetime  cost  of  $1.5  million 


-  [beoff  Andersen,  USAFA 
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Projective  Quantum  Measurements 
on  Spatial  Modes  of  the  Optical  Field 


Quantum  Key  Distribution  (QKD)  for  generation  of 
encryption  keys  over  free-space  optical  links 

Problem: 

•  Optical  losses  in  free-space  quantum  channels  slow  key  generation  rates 

Approach: 

■  Increase  information  carrying  capacity  of  single  photons 

■  Exploit  spatial  attributes  (a.k.a.  linear  and  orbital  angular 
momentum  quantum  states)  of  the  photon  wave  function 

■  Develop  transmission  volume  holograms  as  quantum 
projection  operators  in  order  to  utilize  spatial  modes  in  a 
QKD  protocol. 


VO  •  • 


V 


To  date,  have  demonstrated: 

■  Deterministic  photon  measurements  in  4  dimensions 

■  provides  fidelity  in  QKD 

■  Random  photon  measurements  in  4  dimensions 

■  provides  security  in  QKD 


Experimental  Results: 


spatial  light  modulator 
computer-generated 


Enhancement  to  Bit  Generaton  Rate 


system  Fresnel  number,  D7D2/2Az 
(aperture  sizes  /  propagation  distance) 


volume 

holograms 


single-photon 

detectors 


Mark  Gruneisen,  AFRL/RDSS 


laser 


attenuator 

DISTRIBUTION  STATEMENT  A  -  Unclassified,  Unlimited  Distribution 


Measured  Photon  Detection  Probabilities 


matched  mutually  unbiased 

measurement  basis  measurement  basis 
(deterministic  outcome)  (random  outcome) 


Beam  Control  through  Turbulence 


4i» 


WHftT  is  TueBuuewce. 


fib, 

Wk 
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The  Disturbance 


The  EM  Field 


RDS  -  ASALT 


JTO  -  HTS 


JTOJBCID(tbd) 


AFOSR  MURI 


JTO  MRI 


AFOSR  -  grant 


AEHQR  QTTR 

ACnCR  CTTR 

AFOSR  STTR 

AFOSR  -  RDS  A 


SBIR 


AFOSR  -  RDSM 


SBIR 


AFOSR  -  grant 
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I  Wave  Optics  of  Deep  Atmospheric  Turbulence: 
From  Underlying  Physics  towards  Predictive  Modeling, 

Mitigation 


Co-Investigators  and  Institutions:  M.  Vorontsov  (UD),  S.  Fiorino  (AFIT),  M.  Roggemann  (MTU),  S.  Basu  (NCSU), 
D.  Voelz  (NMSU),  O.  Korotkova  (UM) 


Building  bridges  between  weather  research 
and  forecasting,  computational  fluid 
dynamics  and  statistical  wave-optics 


Atmospheric  optics  effects  sensing  and 
predictive  modeling:  from  Kolmogorov 
idealization  towards  real-world  complexity  | 


Super-resolution  in  volume  turbulence 
Super-focusing  in  deep  turbulence 


ulence-induced 
enhanced  correlation 


Atmospheric  optics 
effects  exploitation: 
considering 
complexity  as  an 
opportunity 


M STATEMENT  A 


•  Sensing  of  large-scale  coherent  atmospheric  structures 

•  Novel  atmospheric  sensing  techniques 

•  Combined  turbulence,  refraction,  &  scattering  effects 


•  Target-in-the-loop  laser-beam  propagation 

•  Propagation  inside  atmospheric  stratified  layers 


Turbulence 

mitigation 


Phase  Fiber 
shifters  amplifiers 


sensing  module 


Sparse  AO 
Exotic  beams 


Controllable 

coherence 


Mitigation  of 
atmospheric  effects: 
foundation  for  intelligent 
engineering  of  optical 
waves  and  systems 


How  Deep  is  "Deep  Turbulence?" 


1 )  Cumulative  effects  of  wave  length,  path  length  and  optical  turbulence  level  (Cn2 ) 

2)  Saturation  of  the  normalized  variance  of  intensity  fluctuations. 

3)  As  intensity  variance  is  saturated,  must  use  phase  information. 

4)  Initially  estimate  both  the  Cn2  and  outer  scale  from  at  short  path  or  using  long 
wavelength  at  low  a  Rytov  value 

5)  Use  this  outer  scale  value  and  the  reconstructed  unwrapped  phases  in  the 
presence  of  branch  cuts  to  estimate  Cn2  at  deep  turbulence  levels  (saturation). 


Physics  of  Distributed  Volume  Turbulence 


U.S.  AIR  FORCE 


©AS ALT 

The  Next  Generation  in  Adaptive  Optics 


Field  tests  corroborate  theory-lab-sim  results 


uib  i  nibu  I  iun  b  i  a  i  bMtN  i  a  —  unciassiTieci,  umimiTea  uistriDution 
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YIP:  SAL  Imaging  and 
Atmospheric  Turbulence 


R  =  yJ(R,  -  y)l  +  (nd  -  x)2 


k-space  representation 


Example  table-top  SAL 
images  made  by  MSU- 
SL  with  broadband 
chirped  lasers 


Atmospheric  Turbulence 


5  10  15  20  25 

Cross-range  [mm] 


Use  high  resolution  FMCW 
ladar  tools  to  study 
synthetic  aperture  ladar 
(SAL)  imaging,  atmospheric 
turbulence,  and  their 
interaction. 


SAL  and  Turbulence  Simulations 


(Left)  Peak  range  vs.  time  shows  several  mm  of 
apparent  movement 

(Right)  FMCW  ladar  history  of  atmospheric 
turbulence  induced  motion  at  14km  slant  range. 


Phase  Screen  for  cj;  =  1e-15  m 


[meters] 


SAL 

Processing 


-0.015  -0.01  -0.005  0  0.005  0.01  0.015 

Cross-Range  [m] 


Added 

Turbulence 


Cn2-1x10‘15  Cn2-1x10'14 


^ElI^ 
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Propagation  of  Polarization-Modulated 
Beams  through  a  Turbulent  Atmosphere 


Amplitude 


ENTERING 

CONTROL  BE  AM. 64'  tun  LIN  POL 


ENTERING 
MODULATED  BEAL  I 
(ft  an  Argon  lme  &  64'  nm 
&  EITHER  POLMOD  or 
AXMOD 


Terrance  Barrett,  BSEI 


POLMOD  achieved  by  combining  two  orthogonally 
polarized  beams  of  different  wavelength. 


VOLTS:  UNMODULATED  CONTROL  BEAM 


“Line  of  No 
Difference” 
Above  this  line 
modulated  beam 
is  higher 
amplitude  than 
control  beam. 


The  control  beam :  647  nm  vertically  polarized  in  this  case. 
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Imaging  through  Turbulent  Media 


Phase  Conjugation  in  Resonantly-Enhanced  Second-Order  Nonlinear  Materials 

Blurred  images  caused  by  atmospheric  turbulence  can  be  completely  restored  to  their  original  resolution 
Measure  reflecting  phase-conjugate  wave,  after  passing  through  turbulent  medium 
Theory:  Y.  J.  Ding,  Opt.  Lett.  37,  4792-4794  (2012) 

Advantages: 

Polarization-insensitive 

Low  pump  powers  required  by  using  second-order  nonlinearity 
Response  is  instantaneous 

-  suitable  for  correction  of  blurred  images  caused  by  atmospheric  turbulence 


Images  before  phase  distortion  for 
(a)  P-polarized  and  (d)  s-polarized 
input  beam 


(c)  and  (f)  corresponding  images 
restored  by  phase-conjugated  beams 


(b)  and  (e)  corresponding  images 
after  phase  distortion 


Yujie  Ding,  Lehigh  University 
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STTRs 


•  Phase  I 

—  Subaperture  Adaptive  Optics  for  directed  energy 
phased  arrays 

—  Complex  field  atmospheric  sensing  suite  for  deep 
turbulence  research 

*  Phase  II 

—  Beam  Control  for  Optical  Phased  Array 
Transceivers 
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i  iNon-resolved  Space  Object  Identification 
Ibex  Software  for  GEO  Characterization 


(  Entry  Point  j 


5 

Parameter  setup 
flrtd  Control 

_ 

12 

Dn  r  .1  mn  +  nr  cnlnri 

and  Control 

Periodicity  and 
Frequency 
Spectrum 
Analysis 
(PAtool) 

Rotation  State 
Analysis 
(ROTtool) 


Initial 

Stability 

Analysis 

(ISAtool) 


Attitude/Shape 

Retreival 

(ASRtool) 


Color 

Legend 

6.1 -level 
algorithm 
funded 
by: 

AFOSR 

(previous) 

AFOSR 

(ongoing) 

DARPA 

(Ibex) 


DISTRIBUTION  STATEMENT  A  -  Unclassified,  Unlimited  Distribution 


DURIP: 

Falcon  Telescope  Network 


Global  Network  of  small  telescopes 

•  Non-resolvable  space  object 
characterization 

•  Satellite  metrics  and  debris  studies 

•  Sensor  networking 
Robotic  and  automated 

13  fixed  and  2  mobile  sites 
Planned  IOC:  2013-2014 


FTN  Observatory  Configuration 


Global  Coverage  (GEO)  Global  Coverage 

(LEO,  1000  km) 

Partners 

•  U.S 


20-inch  Telescope  Truss 


Paramount  ME 


Boltwood  Cloud  Sensor 


•  Colorado  Mesa  University  (Grand  Junction,  CO) 

•  Fort  Lewis  College  (Durango,  CO) 

•  Northeastern  Junior  College  (Sterling,  CO) 

•  Otero  Junior  College  (La  Junta,  CO) 

•  Penn  State  University  (State  College,  PA) 

Chile 

•  Universidad  de  La  Serena  (La  Serena)  & 
Observatorio  Mamalluca  (Vicuna) 

Australia 

•  University  of  Queensland  (Brisbane) 

•  Catholic  Education  Office  of  Western  Australia  and 
Curtin  University  (Perth) 

•  University  of  New  South  Wales  &  EOS  (Canberra) 
Potential:  South  Africa,  Kauai  Community  College, 
hntemationa^gac^niversit^^^^^^^^^^^^^^^ 


& 

i 

12-foot  Clamshell 
Filter  Wheel  Apogee  F47 
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Accurately  Predicting  the 
Location  of  Space  Objects 


Monthly  N  umber  of  Objects  in  Earth  Orbit  by  Object  Type 


Cornerstone  capability  to  all  SSA 


Challenges: 

-  Large  number  of  objects  (>  20,000  of  size 
greater  than  10cm) 

-  Limited  observations  from  multiple 
sensors 

-  Non-Gaussian  errors,  uncertainty 
analysis 

-  High  fidelity  models  representing  space 
weather 

-  Conjunction  analysis  (combinatorial 
problem) 

-  Closely  spaced  objects 

-  Data  association  from  diverse  sensors 

-  Optimal  sensing  for  maximum 
information  gathering  and  hazard 
assessment  in  a  timely  manner 

-  Efficient  computations 
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Accurately  Predicting  the 
Location  of  Space  Objects 


John  Junkins  (TAMU)  -  Insights  to  understand,  quantify, 
and  speed  up  convergence  using  modified 
Chebyshev/Picard  Iteration 


Large  Efficiency 
"Sweet  Spot" 
{relatively  invariant  to 
(interval,  order)}  ^ 


(a)  Computational  Time 


(b)  MCPI  CPU  Time  Efficiency  Sweet  Spot 


/*'***/ 


,0°  oeg'eeN 


Extremely  Large  Accuracy  ^ 
’’V  ]  f  **tSweet  Sjfot”  fj  * 

( for  It  to  1J  significant 
Ai  V  dibits) 


I 

i 


(d)  MCPI  Accuracy  Sweet  Spot 


(c)  MCPI  Significant  Figures 

Figure  2  LEO  Orbit  Computation  Tradeoffs: 


CPU  Time  Versus  Accuracy  Versus  Approximation  Order  Versus  Length  of  Path  Interval 

(Degree  varying  from  40  to  300.  Path  Interval  from  0.01  to  2.2  orbit  periods) 


Aubrv  Poore/Josh  Horwood  (Numeric  Corp.)  -  Highly 
stable,  parallelizable,  super-convergent  family  of  orbital 
propagators  for  orbit  and  uncertainty  propagation 


O  3000 

"CT5 


Q_ 

E  , 


—  Implicit  Runge-Kutta 
q  .  -  -  -  Dormand-Prince  8(7) 

Adams-Bashlorth-Moulton 

•  ■ 

t 

0  - 

■y- 

■  ■ 

0 

0  - 

-  * 

~  .  . . ▲▲ 

■  4  *  " 

Numerica’s  new 
propagator 

0I — 1 - 1 - 1 - 1 - 1 - 1 - 

Dan  Sheeres  (U  Colorado)  -  Correlation  of  and  change 
detection  in  space  objects  using  non-Gaussian  PDFs  for 


hypotheses-free  correlation 
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*  Phase  II 

—  Time  series  prediction  for  satellite  ballistic 
coefficients 

—  Realistic  State  and  Measurement  Error  Uncertainty 
Computation  and  Propagation  for  Space 
Surveillance  and  Reconnaissance 

—  Innovative  Earth  Gravity  Reformulation  and 
Numerical  Integration  for  Responsive  SSA 
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DISTRIBUTION  STATEMENT  A  -  Unclassified,  Unlimited  Distribution 
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DISTRIBUTION  A:  Approved  for  Public  Release 


Questions? 
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